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Abstract 



CN ' We study bound states in a model with scalar nucleons interacting via an 

J^ ■ exchanged scalar meson using the Hamiltonian formalism on the light front. 

In this approach manifest rotational invariance is broken when the Fock space 
is truncated. By considering an effective Hamiltonian that takes into account 
two meson exchanges, we find that this breaking of rotational invariance is 
decreased from that which occurs when only one meson exchange is included. 
The best improvement occurs when the states are weakly bound. 
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I. INTRODUCTION 

Covariant field theory can be coerced into giving solutions to non-perturbative problems 
such as bound states via the Bethe-Salpeter equation. There are technical difficulties in 
arriving at those solutions, but they are manifestly covariant. An alternative approach that 
also provides covariant results is obtained by using a Hamiltonian theory and the correspond- 
ing Schrodinger equation to solve for the bound states. A consequence of this covariance is 
that the commutation relations between the operators that constitute the four-momentum 
(P M ) and the four-angular momentum (M^ u ) give a representation of the Poincare algebra. 
There are many different formulations of Hamiltonian theories, corresponding to different 
forms of dynamics, as discussed by Dirac JIJ. The components of the P^ and M^ u operators 
can be classified as either dynamic (dependent on the interaction) or kinematic. Different 
forms of dynamics result in different separations of P^ and Af" into dynamic operators 
and kinematic operators. The two forms that we will be concerned with are the equal-time 
and light-front formalisms. The equal-time form gives four dynamical operators, while the 
light-front form yields only three dynamical ones. Having a small number of dynamical 
operators is good since this simplifies some calculations. 

Using the canonical coordinate system, with x^ = (x° , x 1 , x 2 , x 3 ) and defining the com- 
mutation relations at equal time (x° = t), we obtain the most familiar Hamiltonian theory, 
an equal-time theory. For this system, rotations and translations are kinematical, while 
the Hamiltonian and generators of boosts are dynamical. The fact that the generators of 
boosts are dynamical implies that when the full Hamiltonian is truncated, it will no longer 
transform appropriately under boosts. This implies that the equal-time Hamiltonian may 
not be well-suited for relativistic problems. 

An alternative coordinate system is obtained by using light-front variables x^ = x° ± x 3 
to construct x^ = (x + , x~,x x , x 2 ). High-energy experiments are naturally described using 
these coordinates, rather than equal-time ones. This is because the front of a beam traveling 
at the speed of light in the (negative) three-direction is defined by a surface where x + is 
constant. The description of the interaction of this beam with a target is much simpler if 
described in terms of light-front variables [0,|3| than if an equal-time variables are used. As 
an example, the Bjorken x variable used to describe deep inelastic scattering experiments is 
simply the ratio of the plus momentum of the struck constituent particle to the total plus 
momentum of the bound state. 

Using the light-front coordinate system and defining the commutation relations at equal 
light-front time (x + = £lf) 5 we obtain a light-front Hamiltonian theory [@,[H|. For this 
choice, rotation about the three axis, translations in the one, two, and plus directions, and 
the generators of boosts in the one, two, and plus directions are kinematical, while the 
Hamiltonian (P~) and rotations about the one and two directions are dynamical. Since the 
boosts are kinematical, the light-front Hamiltonian is useful for relativistic problems, such 
as strongly-bound systems, even when the Hamiltonian is truncated. However, since two 
of the rotations are dynamical, the momentum operator four-vector will not be covariant 
under those rotations when the Hamiltonian is truncated. This lack of manifest rotational 
invariance presents some difficulties for bound-state calculations. 

A key benefit of using the light front is that the vacuum can be very simple. For our 
theory, the masses of the particles are large. All massive particles and anti-particles have 



positive plus momentum, and the total plus momentum is a conserved quantity. Thus, there 
is no condensation, and the vacuum (with p + = 0) is empty, the Fock space vacuum. Thus, 
the Hilbert space for this theory is the Fock space. Diagrams that couple to the vacuum 
vanish, so the number of non-trivial light-front time-ordered diagrams is greatly reduced 
compared to the equal-time theory. 

A. Difficulties with the light front 

An untruncated light-front Hamiltonian will commute with the total relative angular 
momentum operator, since the total momentum commutes with the relative momentum. 
Thus, for the scalar theory we consider here, the eigenstates of the full Hamiltonian will 
also be eigenstates of the angular momentum. However, as mentioned earlier, a Fock-space 
truncation of the light-front Hamiltonian results in the momentum operator four-vector 
losing covariance under rotations. Hence J 2 and the truncated Hamiltonian do not commute 
and this implies that the eigenstates of the truncated Hamiltonian will not be eigenstates of 
the angular momentum. 

How will this violation of rotational invariance affect physical observables? A way to 
observe this violation is to note that on the light front, rotational invariance about the three- 
axis is maintained. This allows us to classify states as eigenstates of J3 with eigenvalues m. 
We compare the energies of states with the same "relative angular momentum" quantum 
number I but different m values. We will define the "relative angular momentum" operator 



in section |HQ . If the Hamiltonian were rotationally invariant, the energies should be the 
same; the breaking of rotational invariance causes the energies to be different ||. 

We expect that the higher Fock-space components of the full Hamiltonian will be small if 
the coupling constant is small enough. Thus, truncation at successively higher orders should 
reduce the violation of rotational invariance of the truncated Hamiltonian. In particular, 
by retaining enough terms in the perturbation expansion of the Hamiltonian, the violation 
of rotational invariance can be reduced to an arbitrarily small amount, provided only that 
a perturbation expansion is valid. Thus, the only real question is: How many terms are 
required? 

B. Outline of the rest of the paper 

The purpose of this paper is to explore how truncation at different orders in the effective 
potential affects the breaking of rotational invariance of the spectra. The model we are 
using is a Lagrangian with neutral scalar particles coupled via a <p 2 x interaction to other 
neutral scalar particles. The Bethe-Salpeter equation for this Lagrangian is discussed, and 
the ladder approximation to the Bethe-Salpeter equation and its spectra are reviewed. We 
also derive the light-front Hamiltonian from that Lagrangian and write out the light-front 
Schrodinger equation for bound states consisting of two particles. An approximation to the 
light-front potential that causes the light-front Schrodinger equation to be equivalent to the 
ladder Bethe-Salpeter equation is discussed. This approximation will be called the uncrossed 
approximation. The one-boson-exchange (OBE) effective potential and two-boson-exchange 



(TBE) effective potential are derived in this approximation. All of this is discussed in more 
detail in section ||. 

Following the development of the model, we look at some numerical calculations in sec- 



tion |T|. For our tests, we look at states composed of two massive scalar particles bound by 
the effects of the exchange of a lighter scalar particle. We calculate the spectra (the coupling 
constant versus bound-state mass curves) using the Bethe-Salpeter equation in the ladder 
approximation. Since the Bethe-Salpeter equation is manifestly covariant, its solutions are 
states with definite angular momentum. We also calculate the spectra using a light-front 
Schrodinger equation with the OBE and the OBE+TBE potentials. Because of the lack of 
manifest rotational invariance, these wavefunctions do not have definite angular momentum. 
An artificial "relative angular momentum" operator is constructed, and a partial-wave de- 
composition is performed on these wavefunctions to analyze the various angular momentum 
states present in them. The problem of trying to classify these states is discussed. We then 
plot the spectra for the wavefunctions classified as the lowest-lying p- and <i-wave states for 
the Bethe-Salpeter, OBE Schrodinger, and the OBE+TBE Schrodinger equations. 

In section |TVJ we summarize our findings: the higher-order terms help to partially restore 
the breaking of rotational invariance introduced by using only the OBE potential. This 
restoration is largest when the states are weakly bound. 

We stress that the point of this paper is not to solve the Lagrangian of this model 
exactly. Instead, we seek to understand how rotational invariance is broken by the light-front 
Hamiltonian in the uncrossed approximation, how it can be partially restored by keeping 
higher-order graphs, and also to compare the results from the two truncated Hamiltonians 
to the results of the "exact" theory (the ladder Bethe-Salpeter equation). Discussion of the 
inclusion of the crossed-graph contributions will be given in a forthcoming paper [|7| . 

This study is closely connected to the work of Schoonderwoerd, Bakker, and Karmanov 
||, who considered two-body scattering in perturbation theory using a model similar to 
ours. Off-shell scattering amplitudes were computed at order g A and it was found that for 
low momenta, the TBE diagrams contributed much less than the iterated OBE diagrams. 
This lead to a conjecture that the lightly-bound states should be well described by the 
OBE potential, with the TBE potential playing a minor role. These results encourage us to 
calculate and compare bound states incorporating the OBE and OBE+TBE diagrams. We 
note that some of the difficulties present in H are avoided here since there are no singularities 
in the bound-state integral equation. 

While this manuscript was in preparation, a paper by Sales et. al. |J appeared. In that 
work, a model similar to ours was used to consider bound states. Their study included 
the same two truncations that we use, as well as the a comparison to the Bethe-Salpeter 
equation results. However, the investigation of Ref. focused on the ground state, while 
here we are concerned with the rotational invariance of the excited states. 

II. OUR MODEL 

We consider two distinguishable uncharged scalars (f) = (</>i, $2) with mass M (which we 
will refer to as nucleons), that couple to a third, uncharged scalar \ with mass \x (which we 
will refer to as a meson) by a (f) 2 x interaction. This <fr 2 x theory, which can be considered 



a massive extension of the Wick-Cutkosky model [jID| , has been used by Schoonderwoerd, 



Bakker, and Karmanov || on the light front to study scattering states. The Lagrangian is 

£ = \ (d^d^ - M V) + \ frxtrx ~ A* V) + 9 -fx<P 2 . (1) 

This Lagrangian will be used in two formalisms, the Bethe-Salpeter equation and the Hamil- 
tonian approach. 

A. Bethe-Salpeter equation 

The Bethe-Salpeter equation |TT|-|r5| provides a way of describing bound states based 
on Feynman propagators and kernels constructed from covariant quantum field theory, and 
thus is manifestly covariant. The equation for the bound state of nucleons 1 and 2 can be 
written as 

GKip = ip, (2) 

G is the free two-particle propagator, which is the product of two one-particle propagators, 
G — S1S2, ip is the four-dimensional Bethe-Salpeter amplitude, and K is the sum of all 
two-particle irreducible two-to-two Feynman graphs, 

— • — • • 4»» — • — — f — 
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K = ! + n + ! + ! + + " ' • ( 3 ) 

— 4— -rf-k- —4— -t4«- — *— 

There are no nucleon exchange graphs since we treat only the case where the two nucleons 
are distinguishable. 

We consider the ladder approximation to the Bethe-Salpeter equation, 

g 2 GK OB ^ = V, (4) 

which is obtained by replacing K in Eq. (Q) with g 2 K BE, the graph due to one-boson- 
exchange, 
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Kobe = — ! • (5) 

This definition of -Kobe makes it independent of the coupling constant. Making this approx- 
imation leaves the Bethe-Salpeter equation covariant. This implies that the Bethe-Salpeter 
amplitudes i[) have definite angular momentum I, and hence the energies of the states are 
degenerate for different m projections of the same angular momentum. 

This equation can be simplified in the center-of-momentum frame. In that frame, once 
the total energy is defined as P°, the four- momentum of the second particle is given by 
&2 = P^ — ki, and thus the Bethe-Salpeter equation is effectively a one-particle equation. 
The remainder of the discussion in this section will be done in the cm. frame. We then can 
rewrite Eq. (|), taking into account the explicit symmetries of the equation, as 

1 1 2 

#lbse(- e G(E)K OBE (E)il) n j !m = ljj n: l : m, (6) 



where E is an arbitrary energy, if} n ,i,m * s the nth Bethe-Salpeter amplitude with angular 
momentum I and three-projection m, and #lbse(^) * s the coupling constant which yields 
that bound-state Bethe-Salpeter amplitude with E as the bound-state energy. We call 
this g{E) the spectrum of the ladder Bethe-Salpeter equation for the corresponding Bethe- 
Salpeter amplitude. The Greens function G(E) and the OBE kernel Kqbe{E) are functions 
of the energy in the cm. frame and are implicitly effective one-particle operators. 

The comparison we draw is only between the ladder Bethe-Salpeter equation and the 
light-front Hamiltonian that corresponds to that approximation. The exact nature of the 
correspondence is discussed in section |U(J| . Since we are not looking at the solutions to the 
full theory, for our purposes it does not matter that there are sizeable differences between 
the solutions to the full Bethe-Salpeter equation and the ladder approximation when the 
coupling constant is large 



B. Light-front Hamiltonian 



To obtain the light-front Hamiltonian from the Lagrangian in Eq. ([3]), we follow the 
approach used by Miller || and many others (see the review 0) to write the light-front 
Hamiltonian (P~) as the sum of a free, non-interacting part and a term containing the 
interactions. This is accomplished by using the energy-momentum tensor in 



p» = - f dx~d 2 x ± T~ 



[x 



\J - Jb - *Xj 



(7) 



The usual relations determine T +M , with 



T HV = _ g „u C + ^ 



dC 



d{d^(j) r j 



-d\ 



(8) 



in which the degrees of freedom are labeled by <p r . 

It is worthwhile to consider the limit in which the interactions between the fields are 
removed. This will allow us to define the free Hamiltonian P ~ and to display the necessary 
commutation relations. The energy-momentum tensor of the non-interacting fields is defined 
as Tq V '. Use of Eq. (|8|) leads to the result 



rriflU 
J 



fUJ 



<9><9> - — 



d n 



M 



2/2 
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fUJ 



+ dr x vx - -j- [d<rXd a X - Sx 



2,. 2 



(9) 



with 



T" 1 

J- n 



+ M 



2 j. 2 



v ± x • v ± x 



2 2 

H X ■ 



(10) 



The scalar nucleon fields can be expressed in terms of creation and destruction operators: 



(x 



d 2 k ± dk + 6(& 



(27r) 3 / 2 V2^ 



a t (k)e- tk - x + aj t (k)e 



ik-x 



(11) 



M 2 +k 2 . 



where i = 1,2 is a particle index, k ■ x = \{k x + + k + x ) — kj_ ■ x± with k 

and k = (k + , k±). Note that k~ is such that the particles are on the mass shell, which is a 
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consequence of using a Hamiltonian theory. The 9 function restricts k + to positive values. 
Likewise, the scalar meson field is given by 



X(x) 



d 2 k ± dk + 9(k + ) 
(2tt) 3 / 2 v / 2A^ 



a x {k)e~ ik - x + aUk)e ik - x 



(12) 



where k 



n 2 +k 2 . 



-p: — , so that the mesons are also on the mass shell. The non- vanishing 
commutation relations are 



a a (k), al(k') = 6(k± - k' ± )6(k + - k' + ), 



(13) 



where a = 1, 2, % is a particle index. The commutation relations are defined at equal light- 
front time, x + = 0. 

The derivatives appearing in the quantity T + ~ are evaluated and then one sets x + to 
to obtain the result 



Pn 



M 2 + k 2 a 2 + k 2 

"- 'a\(k)ai(k) + a} 2 (k)a 2 {k)j + — ± a x (k)a x (k) 



k+ 



(14) 



with J k = J d k±dk + 9{k + ). Eq. ( |i4|) has the interpretation of an operator that counts the 

_ A ^2 + fc 2 U 2 +k 2 

light-front energy k (which is — -^-^ for the nucleons and fc+ x for the mesons) of all of 
the particles. 

We now consider the interacting part of the Lagrangian, Cj. An analysis similar to that 
for the non-interacting parts yields the interacting part of the light-front Hamiltonian Pf; 



Pf=J2 



8=1,2 



g_M_ 

2 Jk Jk 



' (2ir) 3 / 2 s j2k+k>+(k+ + k'+) 
x | 2a\(k + k')a x (k')ai(k) + a x {k + k')ai(k')ai(k) 



-l-Hermitian conjugate} . 



(15) 



For our theory, the Hilbert space is the Fock space. Eq. (|15| ) is self-adjoint on that Hilbert 
space, since the equation is written in terms of the Fock operators. The total light-front 
Hamiltonian is given by P~ = P ~ + Pf. 



C. Hamiltonian bound-state equations 



We will be studying the bound states of two distinguishable nucleons in the light- 
front analogue of old-fashioned (time-ordered) perturbation theory: light-front time-ordered 
perturbation theory (LFTOPT). Using this perturbation theory for our Hamiltonian, we 
can write out a two-nucleon effective Hamiltonian "eigenvalue" equation (a light-front 
Schrodinger equation) where the potential is expanded in powers of the coupling constant 
g. This is not a true eigenvalue equation since the potential is energy dependent. We write 
this in the form of Eq. (ffi) 



>«T + V( 9 Z FSE (P-), P-)} \r) = \r)P- (16) 

oo 

V(g,P~) = j:g 2t V m (P-), (17) 



i=l 



where P~ is an arbitrary light-front energy, |?/> n ) is the the n th wavefunction, and 5 , lfse(-P) 
is the coupling constant which yields that bound-state wavefunction with P~ as the bound- 
state energy. We call this g(P~) the spectrum of the light-front Schrodinger equation for the 
corresponding wavefunction. Only even powers of g appear in V since every meson emitted 
must be absorbed. Also, since V(2%) is the potential due to the exchange of i mesons, we call 
V( 2 ) = Vobe and V( 4 ) = Vtbe- The potentials V^i) can be calculated from the light-front 
time-ordered diagrams. 

We want to approximate our potential V so that Eq. fll6|) is physically equivalent to 
the ladder Bethe-Salpeter equation. This approximation of V will be called the uncrossed 
approximation. By physically equivalent, we mean that the spectra of the potential V should 
reproduce the spectrum for the states of the Bethe-Salpeter equation, excluding the so-called 



"abnormal" states ||10| . It is well known how to reduce the Bethe-Salpeter equation to a 
physically equivalent Hamiltonian (Schrodinger-type) equation. For an extensive discussion 
of this issue of defining the potential equivalent to a sum of Feynman graphs in the equal- 
time case see, for instance, Klein ||20|| , Phillips and Wallace |2l], Lahiff and Afnan |22j . 



and for examples on the light front, Ligterink and Bakker p^j . The general procedure to 
get the effective potential due to n boson exchange takes two steps. First, write the sum 
of all Feynman graphs obtained from iteration of the Bethe-Salpeter kernel with n boson 
exchanges. Then, write that sum in terms of LFTOPT graphs, and discard all graphs which 
are not two-particle-irreducible with respect to the light-front two-body propagator, 

°» ( p ~) = p-hj?- (18) 

The graphs which remain after this procedure constitute the Kbe- 

As an example, we construct the TBE potential. When the ladder Bethe-Salpeter equa- 
tion is used, only one Feynman graph contributes, the box diagram arising from the iteration 
of the Feynman OBE kernel. This gives six non-vanishing LFTOPT diagrams. (The other 
diagrams vanish because the vacuum is simple on the light front.) 



• • 



1 ' ) = / / + / \ + \ / + \ \ + // + \\ . (19) 

4-4-/Feynman -±-4- -4— 4- -44— -^-4- -44 — — W- 

The first four diagrams are iterations of the OBE potential and are reducible with respect 
to Glf- The last two are two-particle-irreducible and thus constitute the TBE potential, 

Vtbe- 

A truncation must be made for the potential in Eq. (|TTD, since in this Hamiltonian theory 
the infinite sum of graphs (required to reproduce the covariant results of the Bethe-Salpeter 
equation) cannot be calculated. In this paper we consider two truncations, one where we only 
keep the OBE potential, and another where we keep OBE and TBE potentials. We write 
out the matrix elements of Vobe and Vtbe in the two-particle momentum basis, denoting 
the momentum of the incoming particles by p\ and P2, and the outgoing particles k\ and &2- 



For Vtbe; particles 1 and 2 have intermediate momenta which we denote by q\ and q 2 . In 
terms of diagrams, the potentials are 

fa m pi fa m p\ 

Vobe = /Pi -fa + \fa~Pi (20) 





/ Pi - fa + 
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P2 
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ft - 


ft pi 

— T^ 

fa // Pi-qi + 


A-2 


q 2 P2 



— •- 


\fa- 


-pi 


fa 


ft 


P2 


fa- 


ft\ 


\ft 



Pi 



^tbe = ft - fci / / Pi - ft + fci - gi \\ ft - Pi • (21) 

WV- 

fc 2 ft P 2 

These potentials are implicitly functions of the incoming and outgoing momenta. 

The wavefunction in Eq. (|16D can be written as a function of the momenta of the two 
particles. However, since the total P± and P + momenta commute with the interaction, 
the total P± and P + of the bound state must be equal to the sum of the two particles' 
perpendicular and plus momentum both before and after the interaction. Therefore, we 
can take the total momentum of the bound state as an external parameter, and solve the 
light-front Schrodinger equation for that total momentum. Thus, the wavefunction, when 
parameterized by the total momentum,is a function of only one momentum variable. We 
can take that momentum to be the momentum of particle 1. The components of particle 
2's momentum are then k 2 = P + — kf and fc 2 _|_ = Pj_ — fei,j_; the minus component 
(the light-front energy) is defined by the requirement that the particle is on mass shell, so 
fej = (M 2 + kl ± )/k+. 

We also define x = kf / P + = xbj, where xbj is the Bjorken x variable, so that k 2 /P + = 
1 — x. Likewise, we write the Bjorken variables that correspond to the other momenta in the 
diagrams in Eqs. ( p0|) and (|2T|) as y = pf / P + and z = ql / P + . Also, a meson with arbitrary 
momentum q + and q± (not to be confused with the loop momentum qi), has light-front 
energy u"{q) = (/i 2 + q\)/q + . 

Inspecting the rules for converting a light-front time-ordered diagram into a formula [0 , 
we find that each piece of the potential is proportional to 

1 (22) 



2(2Tr)3y/x(l -x)y{l-y) 

We will suppress this term from the following potentials. 

For further simplification, we say that we are in the bound state's center-of-momentum 
frame so that the total four-momentum P M = (P~,P + ,Pj_) = (E,E,0±), where E is the 
bound-state energy. We can do this without loss of generality since the boosts are kinematic, 
so we can easily boost the light-front Schrodinger equation to a frame where the bound state 
has arbitrary momentum. All further calculations will be done in the cm. frame. 

Now we use the light-front time-ordered perturbation theory rules to calculate the po- 
tential due to OBE, 



VoBE{ki,±,x,pt !± ,y; E) = (—) 



B(x-y)/\x-y\ 



E -Pi -k 2 -ix) (fa -pi) 



0(y-x)/\y-x\ 



E -k x -p 2 -lo (pi - hi) 



(23) 



We also calculate the potential due to the TBE, 



dz 



(2vr)3 [J, *(l--z)(z-y)(s-*) 

1 



x 



X 



E-q 1 -k 2 -u (k x - q x ) 

1 



E-p 1 -k 2 -u (fci-?i)-w (gi-Pi) 

l 



x- 



■E-.Pl ~?2 - w (?l "Pi; 

+ {1^2}. 



(24) 



The {1 <^ 2} means to replace all labels 1 with 2 and vice versa, as well as replacing the 
Bjorken variables x, y, and z with 1 — x, 1—y, and 1 — z. The angular part the d 2 q± integral 
in Eq. ( p4j) can be done analytically. The radial part of that integral and the z integral need 
to be done numerically. A detailed discussion of the evaluation of Eq. (|24f) will be given in 
Ref. [0. 

We will find it useful to convert from our light-front coordinates (kf, k±) to equal-time 
coordinates k = (k±,k 3 ) by using an implicit definition of k 3 [23 



E(k) 



E 



2E(k) 



E(k) + k 3 



k 2 



(25) 
(26) 



We will refer to the equal-time vector k as the relative momentum of the two-particle system. 
It is worth emphasizing that we call k an equal-time vector, however no simple change of 
variables can produce an equal-time wavefunction from a light-front wavefunction. This is 
simply a useful change of variables. The physics of the light-front remains from the definition 
of the Hamiltonian (P~) and the commutation relations. 

For scattering states, the OBE potential computed for on-shell nucleons and used in the 
Weinberg integral equation [^SJ (which is essentially the scattering analogue of Eq. (TO)) leads 
to manifestly rotationally invariant results when written in terms of the relative momenta 
||. The similarity between that rotationally invariant result and the usual equal-time result 
implies that this equal-time momentum k can be interpreted as the relative momentum of 
the two particles. For bound states, this exact simplification does not occur, since for bound 
states the potential is, of necessity, evaluated off the energy shell (but on the mass shell). 
However, we expect that the OBE potential, written in terms of the relative momentum, 
is approximately spherically symmetric for lightly-bound states. Thus, the wavefunctions 
are approximate eigenfunctions of the "relative angular momentum" , where we define the 
"relative angular momentum" by the operator L = xx k. Our "relative angular momentum 
operator" is not the same as the true orbital angular momentum operator which is obtained 
from the Lagrangian via the energy-momentum tensor in a way similar to the Hamiltonian. 



10 



Now consider the exchange of the particle labels 1 and 2. This causes 

fci,± -► fe 2 ,± = -fei,i (27) 

kf ^k+ = E- kf, (28) 

which means that k 3 as defined in Eq. (|25|) transforms as k 3 — > — k 3 so k — > — fc. Conse- 
quently, exchange of particle labels 1 and 2 is the same as parity. 

Since the two nucleons are identical except for the particle label, the effective potential 
commutes with parity to all orders in g 2 . Furthermore, the light-front Hamiltonian is explic- 
itly invariant under rotations about the three-axis. These considerations allow us to classify 
the wavefunctions as having eigenvalues p of parity (V) and m of the three-component of 
the angular momentum operator (J3). We label the wavefunction as {ki\ip m ^), where 

<*i|.JslC> = <*ilC>m, (29) 

(h\v\r m , P ) = (h\r m , P ) (30) 

= (fcilC> (31) 

With this, we may write the OBE truncation of full uncrossed Hamiltonian as 

P OBE (g, E) = P - + g 2 VoBE(E), (32) 

which gives OBE light-front Schrodinger equation 

PSb^obKE), E) |C'S BE ) = IC'S BE )^, (33) 

where E is an arbitrary energy, |^mp BE ) i s ^ ne nth wavefunction with parity p and J3 
quantum number m, and gQ^(E) is the coupling constant which yields that bound-state 
wavefunction with E as the bound-state energy. 
For the OBE+TBE truncation, we have 

P T - BE (S, E) = P~ + 9 2 Vobe(E) + g 4 V TBE (E), (34) 

which gives TBE light-front Schrodinger equation 

^f BE («fm e) ic T P BE ) = \r m i BE )E, (35) 



where the quantities here are defined in an analogous way to Eq. flSB]). By comparing the 
spectra 5'obe P (-^') an d 9tbe P (E), we can see what effect adding the TBE potential to the 
OBE potential has on the coupling constant for a given bound-state energy. 

D. Comparison 

The solution to the untruncated light-front Schrodinger equation in the uncrossed ap- 
proximation is equivalent to the solution of the ladder Bethe-Salpeter equation. When the 
full uncrossed Hamiltonian is truncated, differences will be introduced. Thus, here we think 
of the ladder Bethe-Salpeter equation as the exact theory which the truncated light-front 
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Schrodinger equations approximate. As more graphs are included in the truncated light- 
front Hamiltonian potential, the agreement with the Bethe-Salpeter equation will obviously 
be better. The question we wish to answer is how well the spectra g(E) for the two different 
truncations [Eqs. ( |33"D and fl35[)] approximate the spectra for the "exact" theory, the BSE 
results. 

The lack of manifest rotational invariance of the truncated Hamiltonian theory causes a 
breaking of the degeneracy of the spectra of the truncated light-front Schrodinger equations 
for different m states — unlike the case for the Bethe-Salpeter equation. The wavefunctions 
from the Hamiltonian approach are classified by their dominant angular momentum contri- 
bution I, so that we can compare the spectra for different m projections of the same total 
angular momentum I. By doing this, we can compare how the degeneracy of the spectra 
is broken for the OBE and the OBE+TBE truncations, and also compare to the spectra 
obtained from the ladder Bethe-Salpeter equation. 



III. RESULTS 

For our numerical work, we pick the meson mass to be 0.15 times that of the nucleon, 
so \x — 0.15M. This is chosen so that our ground state can be considered a toy model of 
deuterium, and also to facilitate comparison with the results of Schoonderwoerd, Bakker, 
and Karmanov ||. 

The technology for doing these bound-state Bethe-Salpeter equation calculations was 
developed over 30 years ago [|l^,|l^,^]. Since the eigenstates of the ladder Bethe-Salpeter 
equation are also eigenstates of the total angular momentum, there is exact degeneracy in 
the energies of the different m states for the same angular momentum. For the range of 
parameters used in this study, we find that the numerical errors in g 2 are less than 0.5%. 
The numerical errors are largest for the most deeply-bound states (E ph 1.85M) with the 
largest coupling constants (g 2 /4:7i ph 50). 




FIG. 1. We display here the first three energy bands for the ladder Bethe-Salpeter equation. E 
is the energy of the bound state of two nucleons, and M is the mass of each of the two nucleons. 
The mass of the meson is fi = 0.15M. 
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The solutions to the ladder Bethe-Salpeter equation form bands where excited states 
with different orbital angular momentum are approximately degenerate with each other, as 
shown in Fig. |I[ This approximate degeneracy is due to the fact that when \x — > 0, this 
model can be shown to have the same degeneracies as the non-relativistic hydrogen atom 
||10|| . Thus, in that limit, all states with the same principal quantum number have the same 
energy. When fi ^ 0, that degeneracy is broken, but only slightly, as we can see from Fig. 
|TJ. Because of this, we will label our states using atomic spectroscopy notation. 

Next, we consider the two light-front Schrodinger equations given by Eqs. (|33|) and 
(|35|). These equations are solved numerically (for each parity and several m values) for the 
spectrum g(E) for a range of bound-state energies E. The symmetries of the light-front 
Hamiltonian allow us to classify the states according to m and the action under parity, so 
as an example, we calculate the spectra with even parity and m = 0. For the range of 
values we use, we find the numerical errors in g 2 are less than 2%. The errors are largest for 
the most deeply bound states with the largest coupling constants, as was the case for the 
Bethe-Salpeter equation. 




FIG. 2. Spectra for the four lowest wavefunctions of even parity and m = (a) for the OBE 
equation (short-dashed line) and the Bethe-Salpeter equation (solid line), (b) for the OBE+TBE 
equation (long-dashed line) and the Bethe-Salpeter equation (solid line). The curves are labeled 
by n, which indicates that the curve belongs to the n eigenvector. The curve for the Is state 
for the Bethe-Salpeter equation and the curve for the first OBE+TBE wavefunction are very close 
together, but distinct. 



We plot both the OBE spectra [in Fig. |(a)] and the OBE+TBE spectra [in Fig. |(b)] 
along with the Bethe-Salpeter spectra for the even parity and m = states. Based on the 
energies, we see that the n = 1 state is expected to be the Is state, the n = 2 state the 2s 
state, and the n = 3 and n = 4 states the 3s or 3d states. We see approximate agreement 
between both of the truncated Hamiltonian results and the Bethe-Salpeter results for lightly- 
bound systems (where E m 2M). 
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Our states are not manifest eigenstates of total angular momentum J 2 . However, if our 
approximation of the Hamiltonian potential was good enough, we would be able to identify 
the n = 3 and n = 4 states of the light-front Hamiltonian calculation with 3s and 3d states of 
the Bethe-Salpeter equation unambiguously. This would determine the angular momentum 
of the light-front Schrodinger equation wavef unctions. From Fig. 0, we see both that the 
addition of the TBE potential brings the n = 3 and n = 4 states closer to the Bethe-Salpeter 
results, and that in the lightly-bound region the n = 3 and n = 4 states can be identified as 
3s and 3<i states. For more deeply-bound states, this identification cannot be made. 

An alternative approach to assigning angular momenta labels would be perform a partial- 
wave decomposition of the wavefunctions in the real angular momentum basis. The wave- 
functions could then be labeled by the angular momentum component which is dominant. 
In light-front dynamics, this is difficult to do since the perpendicular components of the 
real angular momentum operator (J x and J y ) are dynamical, which makes the angular mo- 
mentum operator as complicated as the Hamiltonian. We are encouraged to look for an 
alternative operator which is easier to use, yet approximates the behavior of the real an- 
gular momentum. We shall use the "relative angular momentum", defined by L = x x k, 
where the vectors are equal-time three-vectors. The relative momentum k is defined by 
Eq. ( P5p and x is canonically conjugate to k. The "relative angular momentum" L is used 
to help analyze our solutions, and is not the same as the real angular momentum that can 
be derived from the Lagrangian using the energy-momentum tensor. However, we will see 
that the "relative angular momentum" gives results that are expected from the real angular 
momentum, so our use of the "relative angular momentum" in place of the real one appears 
to be justified. 

A partial-wave decomposition is performed on the wavefunctions (represented in the 
relative momentum basis) to obtain the radial wavefunctions R^ m for all "relative angular 
momentum" states Y^ m . Since the potential is not manifestly rotationally invariant when 
written in terms of the relative momentum, our wavefunctions will have support from many 
different partial waves. Not all partial waves are allowed; only those partial-wave states 
with the same J3 and parity quantum numbers as the wavefunction give non-vanishing 
radial wavefunctions. We have 

00 

l=m 

We define the fraction of the wavefunction with "relative angular momentum" I as 

/*oo 

fP= dkk 2 \RlZ(k)\ . (37) 

The fractions jf are a measure of the amount of "relative angular momentum" state I in 
the n th eigenfunction. 

To illustrate these "relative angular momentum" fractions, we perform this analysis on 
the four lowest coupling constant wavefunctions with even parity and m = for both the 
OBE and OBE+TBE truncations, the same as in Fig. |[ We show /" as a function of E in 
Fig. H (The plots for the m = ±1, ±2 states show similar behavior, with different / values.) 
Note that more higher "relative angular momentum" states contribute for deeper bound 
states. This is the same behavior that would be expected if the real angular momentum was 
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FIG. 3. These are the fractions of each angular momentum state in the four lowest coupling 
constant wavefunctions with m = and even parity plotted as a function of the binding energy 
E/M. The results for the OBE truncation are shown on the left, (a-d), and the OBE+TBE 
truncation are shown on the right, (e-h). Note that the OBE wavefunctions, in general, have more 
support from a larger number of partial waves. 



used to perform the partial wave decomposition. Recall that the real angular momentum 
will commute with the full potential if no truncation is made. However, the potentials we 
use are truncated, and neglecting the higher-order terms breaks the rotational invariance 
of the potential. Also, both the binding and the importance of the higher-order diagrams 
increase with the coupling constant. Thus, both truncations we consider will break rotational 
invariance more when the states are more deeply bound. In the absence of having the real 
angular momentum fractions, we will use the "relative angular momentum" fractions. 

Examination of the fp curves in Fig. ^| shows us that, as postulated earlier, the n — 1 
and n = 2 states shown in Figs. |3](a,e) and |3|(b,f) are predominately s-wave states, which 
we label the Is and 2s states respectively. For the OBE truncation, the n = 3 state shown 
in Fig. |3|(c) is predominately <i-wave (labeled the 3d state) and the n = 4 state shown in 
Fig. |](d) is predominately s-wave (labeled the 3s state), with little mixing between the two 
states. For the OBE+TBE truncation, the n = 3 and n = 4 states are mixtures of both the 
3s and 3d states, as seen in Figs. [|(g,h) In fact, we see a level crossing between the n = 3, 4 
states in that the n = 3 state is predominantly 3d for lightly-bound systems, but as the 
binding increases, the n = 3 state becomes predominantly 3s. Again, this type of mixing 
would also be found if the real angular momentum was used. 
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If this 3s-3d mixing is ignored, then Fig. |] shows clearly that when the TBE is included 
the amount of "relative angular momentum" states mixed in actually decreases as compared 
to the OBE results. This implies that the wavefunctions of the OBE+TBE equation are 
better "relative angular momentum" states than the wavefunctions of the OBE equation. 
Thus the TBE potential restores some rotational invariance to the OBE calculation. 

We attempt to classify the wavefunctions as states with definite angular momentum. 
If a wavefunction has the most support from the "relative angular momentum" state I, 
we classify that state as having angular momentum I. This procedure designates the first 
two wavefunctions for both the OBE and the OBE+TBE cases as s-wave states, which is 
clearly the right thing to do. For the third and fourth wavefunctions of the OBE+TBE 
equation, there are points where the fraction of the s-wave state equals the fraction of the 
d-w&ve state. Here it no longer clear that such a state should be assigned a definite "relative 
angular momentum" value; we do so regardless and analyze the consequences later. 

We now examine the breaking of rotational invariance of the two truncations based on 
the comparison of states with different m projections of the same angular momentum. Since 
the s-wave state has only m = 0, there is nothing to analyse in that case. Further discussion 
of the ground-state s-wave will appear in Ref. JF|. In Figs. f|, ||, and || we plot the Bethe- 
Salpeter bound-state spectra along with the spectra for the states constructed with the OBE 
and the OBE+TBE potentials. The different m states for the Bethe-Salpeter equation are 
exactly degenerate as a result of the rotational invariance of the equation. The curves from 
the Hamiltonian theory do not exhibit the exact degeneracy in m; the degeneracy is broken 
whether the OBE or OBE+TBE potentials are used. However, we see that the degeneracy 
is partially restored when the TBE potential is included, in that for a given binding energy 
the spread of the coupling constants is always smaller. 
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FIG. 4. Here we plot the curves for the first p-wave state, the 2p state. The Bethe-Salpeter (BS) 
result is plotted with a solid line, the OBE results with the short-dashed lines, and the OBE+TBE 
with the long-dashed lines. 



In Figs. £| and |], we plot the spectra for the 2p and 3p states, respectively. In both 
figures, we see that the m = and m = ±1 curves move closer together after addition of the 
TBE potential. However, the average of the two curves for each case does not move much. 
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FIG. 5. Here we plot the curves for the second p-wave state, the 3p state. The Bethe-Salpeter 
(BS) result is plotted with a solid line, the OBE results with the short-dashed lines, and the 
OBE+TBE with the long-dashed lines. 



In fact, for the 3p state in Fig. |5], the m = ±1 curve moves farther away from the ladder 
Bethe-Salpeter curve after addition of the TBE. We also note that in Fig. ^ the spread of the 
spectra is larger than in Fig. |] for both the OBE and the OBE+TBE light-front Schrodinger 
equations. We attribute this to the neglect of the higher-order graphs. Since the states with 
different m values would be degenerate if all the higher-order graphs were included, not 
including them causes a breaking of the degeneracy. Because the coupling constant is larger 
for the 3p states than for the 2p states with the same binding energy, the omission of the 
higher-order graphs causes a larger breaking of the degeneracy for the 3p states. 

We also consider the first d-wave states. We see in Fig. || that, as in Figs. [| and ||, the 
states with different m values move together after addition of TBE. The effect of the level 
crossing of the 3s and 3d states shown in Figs. ||(g) and |](h) is seen here in that there 
are two disjoint curves for m = 0. When the level crossing occurs, at E/M rj 1.96, the 
wavefunction which originally was the 3s state becomes the 3d state and vice versa. 

There are some problems that this level crossing causes. In the intermediate region where 
the 3s and 3d states are about equally mixed, it is probably not physically sensible to call 
the state a 3s or 3d state. However, plotting the curves gives us an indication of what the 
wavefunctions are doing, and for this case we see that the m = curve is always bounded 
by the m = ±1 and m = ±2 curves. A more restrictive classification scheme would leave a 
gap (in E/M) between the two m = curves and it would not be so clear that the bounding 
of the m = curve which we observe here occurs. 

Finally, we note that the spread of the curves in the 3d case is approximately the same 
as the spread in the 3p case for each truncation. This tells us that the higher-order graphs 
have the same qualitative effect in both states. 
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FIG. 6. The first d-wave state, the 3d state. The Bethe-Salpeter (BS) result is plotted with a 
solid line. In (a) the OBE results are the short-dashed lines, and in (b) the OBE+TBE results are 
the long-dashed lines. Note that there are two separate m = curves when the OBE+TBE are 
used, one extending a quarter of the way over, and the other hiding under the m = ±1 curve with 
the OBE+TBE. 



IV. CONCLUSIONS 



In this paper, we considered two truncations of the light-front Hamiltonian, the OBE and 
OBE+TBE truncations. Using these truncations, a "relative angular momentum" operator 
was used to study the partial-wave decompositions of the bound-state wavef unctions. We 
found that the "relative angular momentum" operator acting on those states yield behavior 
similar to that expected from the real angular momentum. This result encourages the use 
of this operator to classify the states according to their angular momentum values I, and 
to study the degeneracy of the spectra for states with the same I value but with different 
projections m. We found less breaking of the degeneracy when the OBE+TBE potential 
was used than when the OBE potential was used. Both of these findings indicate that the 
OBE+TBE truncation of the Hamiltonian breaks rotational invariance less than the OBE 
truncation alone. 

However, there is still some discrepancy between our truncated Hamiltonian spectra 
and the ladder Bethe-Salpeter spectra. In general, our OBE+TBE results for the p- and 
d-wave states show deeper binding than the Bethe-Salpeter results. Not surprisingly, this 
disagreement is largest for the most deeply bound-states where the coupling is largest. This 



difference would be removed if all the higher-order pieces of the potential were included. 
These higher-order pieces are also needed for the full restoration of rotational invariance. 
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